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We calculate elastic scattering of positronium (Ps) by the Xe atom using the recently developed 
pseudopotential method [I. I. Fabrikant and G. F. Gribakin, Phys. Rev. A 90, 052717 (2014)] and 
review general features of Ps scattering from heavier rare-gas atoms: Ar, Kr, and Xe. The total 
scattering cross section is dominated by two contributions: elastic scattering and Ps ionization 
(breakup). To calculate the Ps ionization cross sections we use the binary-encounter method for Ps 
collisions with an atomic target. Our results for the ionization cross section agree well with previous 
calculations carried out in the impulse approximation. Our total Ps-Xe cross section, when plotted 
as a function of the projectile velocity, exhibits similarity with the electron-Xe cross section for the 
collision velocities higher than 0.8 a.u., and agrees very well with the measurements at Ps velocities 
above 0.5 a.u. 

PACS numbers: 34.80.-i, 34.50.-s, 36.10.Dr 


I. INTRODUCTION 

Recently observed similarities between positronium 
(Ps) scattering and electron scattering from a number 
of atoms and molecules [1-3] in the intermediate en¬ 
ergy range were explained [4, 5] by the dominance of 
the electron exchange interaction with the target atom 
or molecule. An explicit proof of this equivalence was 
given using the framework of the impulse approxima¬ 
tion [4], valid above the Ps ionization threshold. How¬ 
ever, at lower energies the impulse approximation breaks 
down and more sophisticated methods are required. The 
close-coupling method, which includes an expansion of 
the total wave function over the states of the target and 
the projectile, is very challenging computationally. So far, 
such calculations have been carried only for simple tar¬ 
gets like the hydrogen and helium atoms, often using only 
a small number of states [6-9] . 

Recently we developed a pseudopotential method [5] in 
which a nonlocal Ps-atom potential is constructed based 
on the electron-atom and positron-atom scattering phase 
shifts. This method was successfully applied to the calcu¬ 
lation of Ps scattering from Ar and Kr, and gave results 
in good agreement with those of the beam experiments 
[!]■ 

In the present paper we complete our theory for heavier 
rare-gas atoms by performing calculations of Ps scatter¬ 
ing from xenon. An interesting aspect of this problem 
is the question of existence of the Ramsauer-Townsend 
(RT) minimum in the scattering cross sections. It is well 
known that the RT minimum does exist in electron scat¬ 
tering from Ar, Kr, and Xe. However, our previous calcu¬ 
lations [5] did not find it in Ps-Ar and Ps-Kr scattering. 
We explained this by the relative weakness of the van der 


Waals interaction between Ps and a neutral atom as com¬ 
pared to the polarization interaction in electron-atom or 
positron-atom scattering. This results in positive scat¬ 
tering lengths for Ps-Ar and Ps-Kr collisions, in contrast 
to the negative scattering lengths in e^-Ar and e^-Kr 
scattering. Moreover, instead of the RT minimum, we 
obtained what could be called the “anti-Ramsauer max¬ 
imum”, due to the fast increase of the S- and P-wave 
contributions to the elastic scattering cross section. 

The present calculations confirm the above observa¬ 
tions for Ps-Xe collisions. We also extend the earlier work 
[5] by developing a method for the calculation of the 
ionization cross sections for Ps collisions with rare-gas 
atoms, based on the binary-encounter approach. The re¬ 
sults for the total cross sections agree very well with mea¬ 
surements above the Ps break-up threshold. However, re¬ 
cent beam measurements [10] of Ps scattering by Ar and 
Xe show that the cross section decreases towards lower 
energies, which is not supported by the present calcula¬ 
tions. 


II. Ps-Xe COLLISIONS AND A SUMMARY FOR 
THE HEAVIER RARE-GAS ATOMS 

The pseudopotential method for the calculation of elas¬ 
tic Ps-atom scattering was developed and described in 
detail in Ref. [5], and only a brief account is given here. 
We first calculate the positron-atom and electron-atom 
scattering phase shifts in the static (for e"*") and static- 
exchange (for e~) approximations. We then construct a 
local positron-atom pseudopotential and an Z-dependent 
electron-atom pseudopotential. The latter contains a re¬ 
pulsive core which allows one to decrease the overlap be- 
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FIG. 1. Electron and positron s-, p-, and d-wave phase shifts 
for scattering from Xe calculated for in the static (positron, 
dot-dashed red lines) and static-exchange (electron, solid blue 
lines) approximations. The phase shifts obtained using the 
electron pseudopotential, Eq. (2), with the parameters given 
in Table I, are shown by dashed blue lines. The pseudopoten¬ 
tial phase shifts for the electron d wave and all positron partial 
waves are indistinguishable from the static-exchange/static 
phase shifts on the scale of the graph. 

tween the wave function of the scattered electron and 
the occupied atomic orbitals [11]. Atomic units are used 
throughout. 

Figure 1 shows the s-, p-, and d-wave phase shifts for 
positron and electron scattering from the ground-state 
Xe atom described in the Hartree-Fock approximation. 
The pseudopotential for the positron is chosen in the form 

V,(r} = (1) 

With the parameters Zp and Op chosen as indicated in 
Table I, it gives phase shifts that are indistinguishable 
from the static positron phase shifts on the scale of the 
plot. 

The electron pseudopotential is chosen as 

= + ( 2 ) 

where the second term in Eq. (2) represents the repulsive 
core and exchange. The parameters Ze, ae, B, n, and /3 
are adjusted to obtain the best fit of the Hartree-Fock 
phase shifts (taken modulo tt) . In the present calculation 
we have chosen n = 2, with all other parameters listed 
in Table I. Note that the electron pseudopotential is I 
dependent. This is required to effectively describe the ef¬ 
fect of the Pauli exclusion principle due to the occupied 
ground-state electron orbitals of the Xe atom. The corre¬ 
sponding phase shifts are shown in Fig. 1 by thin dashed 
lines. 

After averaging the sum of the positron and electron 
pseudopotentials over the electron charge distribution in 


TABLE I. Parameters of the positron-Xe and electron-Xe 
pseudopotentials. 


projectile 

1 

^p,e 


B 

P 

e+ 

0-4 

25.09 

1.568 

- 

- 

e~ 

0 

54.00 

2.8522 

129.790 

1.53130 


1 

54.00 

2.0522 

50.1760 

0.88638 


2 

24.713 

1.0731 

4.3433 

0.54959 


3 

15.163 

1.2381 

-3.6233 

1.0710 


4 

14.792 

1.3071 

-3.7086 

1.0999 


Ps(ls), we obtain a nonlocal central potential which de¬ 
scribes the Ps-Xe interaction in the static approximation. 
We then add the van der Waals interaction in the form 


where Cq is the van der Waals constant and Rc is a cutoff 
radius. The Cq constant is calculated using the London 
formula [12], Cq = 240.6 a.u., which is accurate to 5% 
[13]. In previous calculations the cutoff parameter Rc was 
varied between 2.5 and 3.0 a.u. with insignificant change 
in the results for cross sections. In the present calcula¬ 
tions for Xe we chose i?c = 3 a.u. 

The integro-differential radial equation for the wave 
function of the Ps centre-of-mass motion (Eq. (16) in 
Ref. [5]) is solved iteratively. With a suitable choice of 
the local part of the interaction potential, this process 
converges quickly, and the solutions yield the Ps-atom 
scattering phase shifts. These are shown in Fig. 2 for the 
three lowest Ps partial waves, L = 0, 1, and 2. As in 
the case of Ar and Kr [5] , inclusion of the van der Waals 
attraction gives a positive contribution to the scatter¬ 
ing phase shifts. For the S and P waves this leads to 
a decrease in the scattering cross sections. This effect 
of virtual excitations of the target and projectile (which 
underpins the van der Waals interaction) was discussed 
previously in Ref. [14]. Similarly, the scattering length for 
Ps-Xe scattering obtained with the van der Waals poten¬ 
tial, A = 2.45 a.u., is smaller than the value A = 3.57 a.u. 
obtained in the static approximation. The latter value is 
in a reasonable agreement with A = 3.77 a.u. obtained 
by Blackwood et al. [15, 16] in the static-exchange ap¬ 
proximation. We see that the Ps scattering length for 
Xe is greater than those for Ar and Kr, which confirms 
our prediction [5] of the growth of the positive scatter¬ 
ing length with atomic number Z. Although the van der 
Waals interaction (which makes A smaller) increases with 
Z, the effect of the Pauli repulsion for heavier atoms is 
stronger. Mitroy and Bromley [17] used the stochastic 
variational method with model polarization potentials 
for the electron- and positron-atom interactions, and ob¬ 
tained values of the Ps-Xe scattering length in the range 
1.50-2.60 a.u., with the recommended value of 2.29 a.u., 
in close agreement with our value. 


Vw{R) = 1-exp 
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FIG. 2. Elastic scattering phase shifts for Ps on Xe, obtained 
with the static-field pseudopotential (solid lines), and with 
the van der Waals potential added (dashed lines). 


Turning to higher energies, we observe an unusual in¬ 
crease of the phase shifts at u > 1.3 a.u. The S'-wave 
phase shift also exhibits a small kink between v = 1.5 
and 1.6 a.u. This behaviour reflects the influence of the 
repulsive wall in the pseudopotential, and in fact the well- 
known failure [11] of the pseudopotential approach at 
higher energies. However, the S-wave and H-wave con¬ 
tributions are small in this energy region (^60 eV) and 
do not influence the behaviour of the total cross section. 

The cross sections for Ps-Xe scattering are shown in 
Fig. 3. In the low-velocity range, the elastic cross section 
in the static approximation (i.e., not including the van 
der Waals interaction) is close to the corrected static- 
exchange results of Blackwood et al. [15, 16] (not shown 
on the graph). However, at higher velocities, in the range 
V = 0.4-1.0 a.u., our cross section decreases more rapidly 
and is substantially lower than that of Blackwood et al. 
[15, 16]. 

The present total scattering cross section is shown in 
Fig. 3 by the dashed red line. It was obtained by adding 
the elastic and ionization cross sections, the latter calcu¬ 
lated using the binary-encounter method (see Sec. HI). 
As was shown before [4, 15], elastic scattering and Ps ion¬ 
ization are the two major processes contributing to the 
total scattering cross section. 

In Fig. 3 we also compare the total scattering cross 
section with the measurements and the e“-Xe cross sec¬ 
tion. The agreement with the experimental data [3] is 
remarkably good. Note, however, that recent measure¬ 
ments below the Ps ionization threshold [10] indicate that 
the cross section continues to decrease towards lower ve¬ 
locities, in contrast with our prediction of a maximum in 
this energy region. To analyze the similarity between e~- 
Xe and Ps-Xe scattering, we also present the e“-Xe total 



projectile velocity (a.u.) 


FIG. 3. Gross sections for Ps and electron collisions with Xe 
atoms as functions of the projectile velocity: solid black line, 
present elastic Ps-Xe cross section; dashed-dotted magenta 
line, elastic cross section obtained in the static approxima¬ 
tion (i.e., without the van der Waals interaction); dashed red 
line, total Ps-Xe cross section; dotted blue line, e“-Xe elastic 
scattering cross section taken from the calculations [18] and 
measurements [19]; solid squares, measured Ps-Xe total cross 
section [3]. 


cross section calculated by Sin Fai Lam [18] for E < 30 
eV {v < 1.485 a.u.). At higher velocities we show the 
measured cross section from Ref. [19], which agrees very 
well with the calculations [18] below 30 eV. The Ps-Xe 
cross section remains substantially lower than the corre¬ 
sponding e“-Xe cross section for velocities up to 1 a.u. 
This makes the Xe case somewhat different from those 
of Ar and Kr, where the proximity of the electron and 
Ps scattering cross sections was observed right from the 
ionization threshold v = 0.5 a.u. 

Figure 4 presents a comparison of the calculated Ps to¬ 
tal cross sections for all heavier rare-gas atoms. They ex¬ 
hibit the same features: the “anti-Ramsauer” maximum 
at low velocities and a broader maximum in the region 

ft! 1 a.u. It is interesting that in a weaker form this 
feature is also present in the close-coupling Ps-hydrogen 
cross section when excitations of both Ps and the target 
(which account for the van der Waals interaction) are in¬ 
cluded (Fig. 5 in Ref. [8]). Compared with Ar and Kr, 
the absolute magnitude of the Xe cross section is sub¬ 
stantially higher. Note also that the Ar cross section is 
initially slightly decreasing, indicating a relatively weaker 
role of the Pauli repulsion in this case. However, due to 
the positive sign of the scattering length, the shallow min¬ 
imum at u = 0.085 a.u. is not of the same origin as the 
true RT minimum. 

In Fig. 5 we present the momentum-transfer cross sec¬ 
tions for Ar in the low-energy region relevant to experi¬ 
ments [20, 21] on Ps thermalization in Ar. Although the 
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Ps velocity (a.u.) 


FIG. 4. Calculated Ps total scattering cross sections for the 
heavier rare-gas atoms. Data for Ar and Kr are taken from 
Ref. [5]. 



FIG. 5. Momentum-transfer cross section for Ps collisions with 
Ar at low energies: solid line is the calculation from Ref. [5]; 
solid squares with error bars are the experimental data from 
Ref. [20] {E = 0.3 eV) and Ref. [21] {E = 1 eV). 


experimental error bars are large, the results of measure¬ 
ments are not inconsistent with the the tendency of the 
cross section to decrease in the energy range between 0 
and 1 eV. Our cross sections are also compatible with the 
measurements of Coleman et al. [22], who obtained the 
mean value 6.4x 10“^® cm^ in the energy range between 0 
and 6.8 eV, although it seems that our values are some¬ 
what too high in this energy range. The fact that the 
cross sections are likely smaller in this region is also con¬ 
firmed by the recent beam measurements of Brawley et 
al. [10]. 


III. BINARY-ENCOUNTER APPROXIMATION 
FOR Ps IONIZATION 

Ps ionization is one of the two main scattering channels 
in Ps collisions with neutral targets. At low collision en¬ 
ergies comparable to the ionization threshold, the elastic 
scattering dominates, but for energies of about 100 eV 
the ionization cross section becomes comparable to or 
even greater than the elastic cross section. Calculations 
show that Ps ionization is the major inelastic channel in 
Ps collisions [4, 15]. Hence, the sum of the elastic and Ps 
ionization cross sections provides a good approximation 
for the total cross section. 

Owing to the diffuse nature of the Ps atom, the cross 
section of its ionization (or break-up) in collisions with 
atoms and molecules can be calculated using the impulse 
approximation [23, 24]. There are two difficulties associ¬ 
ated with such calculations. One is related to the ambigu¬ 
ity of the on-shell reduction of the electron and positron 
scattering amplitudes. The electron and positron scatter¬ 
ing amplitudes which contribute to the total amplitude 
in the impulse approximation are “off the energy shell”, 
since they need to be evaluated for unequal initial (vi) 
and final (u/) electron or positron velocities, and energy 
E ^ vf j 12. The on-shell reduction of Starrett et al. [24], 
which was also used in Ref. [4], assumes that the ampli¬ 
tude is a function of the effective velocity v = max(ui, Vf) 
and the momentum transfer q is linked to the electron 
(positron) scattering angle dsc by g = 2usin(0sc/2). In 
principle, one can consider different methods of on-shell 
reduction, leading to different results for the Ps scatter¬ 
ing amplitude and cross section. The other difficulty is 
computational. In order to avoid lengthy calculations, 
Starrett et al. [24] used the so-called peaking approxi¬ 
mation, which assumes that the Ps wave function in mo¬ 
mentum space varies much faster than the scattering am¬ 
plitude. While this assumption is justified, it introduces 
additional uncertainty in the results. 

In this section we show that one can use a similar but 
simpler approach based on the binary-encounter approx¬ 
imation [25, 26] and avoid both difficulties in the calcula¬ 
tions of Ps ionization. This approach employs the differ¬ 
ential cross sections rather than amplitudes, for electron 
and positron scattering from the target atom. 

Consider the process of Ps break-up in collision with a 
neutral target B: 

B + Ps ^ B + e'^ + e~. 

Assuming that at the instant of collision the electron and 
positron inside the Ps atom are quasi-free, the ionization 
rate due to e~-B collisions is [25] 

Pion = Nb(\v -Vb\ [ da\ (4) 

\ Jae>i / 

where Nb is the number density of particles B, is the 
relative collision velocity, v is the electron velocity rela¬ 
tive to the Ps center of mass, da is the differential cross 
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section for e~-B elastic scattering, and the integration is 
restricted to the angles which result in the energy transfer 
to the electron AE greater than the Ps ionization poten¬ 
tial / = 6.8 eV. The averaging denoted by (• • •) is over 
the electron velocity distribution in Ps. A similar expres¬ 
sion can be written for the contribution, and the total 
ionization rate is found by adding the two contributions 
and neglecting the interference between them. 

Dividing Eq. (4) by the flux density of incident parti¬ 
cles B, we obtain the total ionization cross section due 
to electron interaction with B as 

|v-vb| f da\. (5) 

Jae>i / 

In the laboratory reference frame the heavy particle B is 
at rest, and as a result of scattering, the electron velocity 
changes from u = v — vb to u', |u'| = |u|. The change 
of the electron kinetic energy in the Ps frame then is 

AE = i [|u' -1- VBp - |u -I- VBp] = vb • (u' - u). 

If we direct vb along the z axis and introduce spherical 
angles {0,cj)) and for the vectors u and u', we 

obtain 



AE = vbu{cos 9'— cos9). (6) 

For the ionization process, integration over 9' is subject 
to the restriction 


I <AE< v%, (7) 

where the upper limit follows from the Ps kinetic energy 
in the laboratory frame, consistent with the threshold for 
the ionization process, Mv%/2 > I, M = 2 being the Ps 
mass. With the help of Eq. (6), the constraints (7) define 
the region in the {9, 9') plane: 

cos 9 -\—— < cos 9' < cos 0 -I- —. (8) 

VbU u 

The electron differential scattering cross section from 
a spherically symmetric target B is 

^ = Y,{21 + l){2l' + l)ft,fiPi>{cos9s)Pi{cos9s), 
w 



FIG. 6. Ionization cross section for Ps collisions with heavier 
rare-gas atoms calculated in the binary-encounter approxi¬ 
mation using the electron and positron phase shifts from the 
polarized-orbital method [27-30]. 


This is a five-dimensional integral with respect to the 
variables 9, (p, 9', p', and u. Using the addition theorem 
for the spherical harmonics and writing 


Yimiu) 


Qimicos 9)—^, 

V zTT 


where 0;m(cos 9) are the normalized associated Legendre 
functions, we can perform integration over the azimuthal 
angles (p and p' with the result 


47r p , 3 

Tion = - / duu 

VB Ji/2vb 

pi —I jvBU 


X J d{cos9)\gis{u^ + v% + 2uvb cos9)\ 

X X! /i*(w)/i(w)0i'm(cos0)0/„(cos6') 


^COS 0-\-Vb l'^ 


/ COS 9-\-I/vb 1 


d(cos0')0;m(cos 6»')0i'm(cos 9'). (9) 


where 9^ is the scattering angle in the laboratory frame, 
i.e., the angle between u and u', and // is the scattering 
amplitude for partial wave Z, 

1 _ pii&i(u) 

defined by the phase shift 5i{u). According to Eq. (5), 
the differential cross section should be multiplied by 
|v — vb| = u, integrated over the scattering angles, and 
averaged over the electron velocity distribution in the 
ground-state Ps, 

J_| 2^|2 _ 1 256 

^ 47r 7r(4u^-|-1)^ 


The integration limits above follow from the restrictions 
(8). The positron contribution has a similar form, with 
fi{u) being the positron scattering amplitudes from B. 

To determine the Ps ionization cross sections for Ar, 
Kr, and Xe, we employed the electron and positron phase 
shifts calculated by McEachran and Stauffer using the 
polarized-orbital approximation [27-30]. Figure 6 shows 
the ionization cross sections for Ps collisions with Ar, 
Kr, and Xe over the range of velocities corresponding to 
energies from threshold to 435 eV. All the cross sections 
peak at a relatively low velocity, slightly above 1 a.u., 
and then decrease rather slowly for Ar and Kr, while for 
Xe the cross section drops and then has a broad second 
maximum at u « 3 a.u. 
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Figure 7 compares the present Ps ionization cross sec¬ 
tions for Xe, Ar, and Kr with those obtained by Starrett 
el al. [24] using the impulse approximation. The agree¬ 
ment is generally good, especially for Ar. Note that the 
impulse approximation cross sections for Ar and Kr both 
peak at approximately 7.5 x 10“^® cm^, while for Xe it 
peaks at approximately 10 x 10“^® cm^. On the other 
hand, the binary-encounter cross sections show a progres¬ 
sive increase of the maximum for Ar, Kr, and Xe, which 
appears to be physically reasonable: for heavier atoms the 
electron elastic scattering cross sections are higher, which 
should lead to higher Ps ionization cross sections. It is in¬ 
teresting that for Xe both the present binary-encounter 
approximation and the impulse approximation cross sec¬ 
tions have the second maximum. Its position and magni¬ 
tude differ slightly between the two calculations. 

There might be several reasons for the discrepancies 
observed in Fig. 7. 

1. The electron and positron scattering amplitudes 
used by Starrett et al. [24] were calculated 
used the static-exchange approximation, while the 
polarized-orbital phase shifts of McEachran and 
Stauffer used in the present work take target po¬ 
larization into account. 

2. The impulse approximation calculations of Starrett 
et al. [24] uses on-shell reduction of the scattering 
amplitude, which is not a unique procedure. At the 
same time, the impulse approximation takes into 
account the interference between the electron and 
positron contributions, which is neglected by the 
binary-encounter method used in the present work. 

3. Starrett et al. [24] used the peaking approximation, 
which neglects the velocity dependence of the scat¬ 
tering amplitude on the scale of the velocity spread 
of the Ps wave function in momentum space. This 
approximation might become less reliable at higher 
energies. On the other hand, the impulse approxi¬ 
mation takes into account the Coulomb interaction 
within the electron-positron pair in the final state, 
while the binary-encounter neglects it. 

In view of all these different approximations made in the 
two methods, the agreement observed in Fig. 7 looks very 
satisfactory. 

Experimental data for Ps collisions with Ar, Kr, and 
Xe atoms [3] do not indicate a second maximum or 
plateau in the total cross section as a function of Ps 
velocity. However, the measurements do not go above 
V = 2 a.u., i.e., they perhaps have not reached the regime 
where the ionization cross section dominates the total. 
Note also that the velocity dependence of the measured 
total Ps-He and PS-H 2 cross sections [31] becomes quite 
flat for velocities between 1.5 and 2 a.u. 



Ps velocity (a.u.) 



Ps velocity (a.u.) 

FIG. 7. Ionization cross section for Ps collisions with (a) Xe, 
(b) Ar, and (c) Kr. Black solid line, present calculations; red 
dashed line, impulse approximation calculations of Starrett et 
al. [24]; blue solid line, the contribution to the ionization 
cross section (shown for Xe only). 


IV. CONCLUSIONS 

Interaction of Ps with atoms is mostly controlled by the 
exchange interaction between the electron in Ps and the 
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target electrons, and by the van der Waals interaction. 
For collision energies above the Ps ionization threshold 
the exchange interaction dominates, making the Ps-atom 
cross section look like e“-atom cross section when plotted 
as a function of the projectile velocity. These features are 
described very well by the present approach which com¬ 
bines the pseudopotential method for elastic Ps scatter¬ 
ing and the binary-encounter approximation for Ps ion¬ 
ization (break-up). New calculations for Xe confirm the 
experimental observations [1-3] of similarity between Ps- 
atom and electron-atom scattering. On the other hand, 
recent measurements [10] for Ar and Xe at low velocities 
do not confirm our predictions of the low-energy peak in 


the Ps-rare-gas-atom cross sections. It is possible that the 
pseudopotential model overestimates the P-wave contri¬ 
bution at low energies, and more theoretical work is nec¬ 
essary to describe Ps scattering in this region accurately. 
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